arXiv:1508.01021vl [hep-lat] 5 Aug 2015 


PROCEEDINGS 

OF SCIENCE 



Nucleon-pion-state contributions in the 
determination of the nucieon axiai charge 


Oliver Bar* 

Fakultdt fur Physik 
Universitdt Bielefeld 
Universitdtsstrafe 25 
D-33615 Bielefeld 
Germany 

E-mail: obaer@physik . hu-berlin . de 

The nucleon-pion-state contributions to QCD 2- and 3-point functions used in the calculation of 
the nucleon axial charge are studied in chiral perturbation theory. For sufficiently small quark 
masses and large volumes the nucleon-pion states are expected to have smaller total energy than 
the single-particle excited states. To leading order in chiral perturbation theory the results do 
not depend on low-energy constants associated with the interpolating nucleon fields and apply to 
local as well as smeared interpolators. The nucleon-pion-state contribution is found to be at the 
few percent level. 


The 33rd International Symposium on Lattice Field Theory 
14 -18 July 2015 

Kobe International Conference Center, Kobe, Japan* 


* Speaker. 


(c) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








Nucleon-pion-state contributions to the nucleon axial charge 


Oliver Bar 


1. Introduction 

In recent years quite a few collaborations started to pursue lattice QCD simulations with pion 
masses close to or at the physical value. For example, at this conference the PACS collaboration 
reported their latest results obtained on a 96^ lattice with Mj^ ss 147MeV and Mj^L ss 6 [1]. Such 
simulations require no or only a short chiral extrapolation, so the uncertainties associated with 
this step are essentially eliminated. On the other hand, some complications get more severe the 
smaller the pion mass is. For instance, the signal-to-noise problem [2] gets worse and prevents 
large euclidean time separations in many correlation functions. In addition, the smaller the pion 
mass the more pronounced is the contamination due to multi-particle-states in correlation functions. 

To discuss this second aspect consider the correlation functions needed to compute the nu¬ 
cleon axial charge gA- One usually computes the ratio = C^{t,t')/C 2 {t) of the 3- and 2-pt 

functions 

03(^4') = , ( 1 . 1 ) 

C 2 {t) = Idhrap{Npix,t)Na{0,0)) . (1.2) 

Na,Np are nucleon interpolating fields placed at source time 0 and sink time t. A\ denotes the 
third isospin component of the axial vector current at insertion time t', and r',r are appropriately 
chosen Dirac matrices (see Ref. [3], for example). Using the spectral decomposition one finds for 
t > t' > 0 

^ gA + bxe-^'^^-'''> +bie-^'‘' +cxe-^'^ + .... (1.3) 

The deviations from a constant are due to excited states with the same quantum numbers as the 
nucleon. The exponential suppression is governed by the energy gap AEi = E\ —M^ where E\ is 
the energy of the first excited state. The nature of this state depends strongly on the pion mass and 
the extent L of the finite spatial box realized in the simulation. For pion masses close to the physical 
value and for sufficiently large L one expects the nucleon-pion state with back-to-back momentum 
\p\ = Itz/L to be the one with the smallest energy, Ei En + Eji. Even the three particle state 
with a nucleon and two pions at rest has a smaller energy than the first single-particle excited state. 
Assuming the aforementioned values ss 147MeV and Mj^L ss 6 one expects five nucleon-pion 
states with an energy between 4% and 1.5Ma- 

The size of the excited state contributions does depend on the size of the coefficients b\,b\,ci 
too. These coefficients are ratios of matrix elements involving the excited states. As pointed out in 
Ref. [4] these coefficients can be computed in chiral perturbation theory (ChPT). This idea is not 
new and has been independently put forward in [5, 6, 7] using heavy baryon ChPT. Here we report 
our results using the covariant ChPT formulation, taking into account the chiral expansion of the 
nucleon interpolating fields. Moreover, we compare the size of the nucleon-pion-state contributions 
in three different methods to compute the axial charge. 

2. Setup: Covariant Baryon ChPT 

The framework for our calculation is the covariant formulation of Baryon ChPT [8, 9] (see 
also the monograph [10]). Both the chiral effective Lagrangian and the expressions for the vector 
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Figure 1: Feynman diagrams for the 3-pt function. Solid and dotted lines correspond to nucleon and pion 
propagators, respectively. Squares, diamonds and circles stem from the interpolating nucleon fields, the 
current and the interaction vertex. 


and axial vector currents can be taken from the literature. The chiral expressions for the nucleon 
interpolating fields are less known but were derived in [11, 12]. Starting point are the familiar 
local nucleon interpolating fields in QCD, Ni = {qq)q and N 2 = (^ 75 ^) 75 ^, wifh q denoting fhe 
quark double! confaining fhe up and down quark, and q = q^Cjs (/CJ 2 ). Based on fhe fransformation 
properfies under chiral symmetry and parity [13] these interpolating fields are mapped fo ChPT 
by writing down the most general expression in B ary on ChPT that has the same transformation 
properties. The resulting expression is then expanded in powers of the pion fields as usual. 

Local inferpolating fields are rarely used in lattice simulations. Instead, smeared interpolating 
fields made from smeared quark fields are usually employed, where qsm{x) = / d^yK{x—y)q{y) 
with some gauge covariant kernel K. The details of the kernel depend on the paiticular smeaiing 
(e.g. Gaussian smeaiing [14, 15], the gradient flow [16], efc.) For what matters here only two 
properties are important: i) the kernel is diagonal in spinor space and ii) it is essentially zero for 
distances |x —y| larger than /^smeai-j a “smearing radius”. In that case the smeared interpolating 
nucleon fields have the same transformation properties as their unsmeared counterparts. If, in ad¬ 
dition, the smearing radius is much smaller than the Compton wave length of the pion the smeared 
inteipolating nucleon fields are mapped onto the same pointlike expressions in ChPT as the local 
ones, the only difference is the different low-energy constants (LECs) in these expressions. 

The calculation of the correlation functions in ChPT is straightforward. To LO it involves 
sixteen diagrams for the 3-pt function, depicted in fig. 1. Four diagrams confribufe fo fhe 2-pt 
function, which was already computed in [12]. A finite box with spatial extent L and periodic 
boundary conditions was assumed. The time extent was taken to be infinite, for simplicity. In 
addition, we assumed isospin symmetry. 

3. Nucleon-pion contribution to 

Taking into account the nucleon-pion state contribution the result for the ratio R{t,t') assumes 
the general form 

^ + +bne-^"^' . (3.1) 
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Figure 2: The coefficients and c^'' as a function of and Mt^L for the lowest two momentum 

states (n = 1 and n = 2). 


The sum runs over the discrete momenta allowed hy the periodic boundary conditions. The energy 
gaps are given hy AEn = +Eji{pn) —Mj^. To leading order (LO), the order we are working 

here, we find bn = bn- This is expected if the same interpolating field is used at source and sink. 
Deviations from bn = bn show up at higher order in the chiral expansion. 

The coefficients b„,bn,Cn are dimensionless and can he written as functions of the following 
dimensionless parameters: /L, gA, Mji/Ma and MjiL. Two LO LECs enter, / and g\, hut the 
coefficients do not depend on the LECs associated with the nucleon interpolating fields. The reason 
is chiral symmetry. It dictates that these LECs are overall factors in both the 3-pt and 2-pt function, 
thus they cancel in the ratio R. Therefore, LO ChPT makes a rather definite prediction for the ratio. 
Eor our purposes we can approximate / and gA by the experimentally well-known measured values 
for the pion decay constant and the axial charge. Once this is done the coefficients depend on two 
unknowns only, Mji/Ma and MjiL. 

It is useful to separate a trivial prefactor and the non-trivial ChPT result in the coefficients by 
writing 


bn — 




r3pt 


m„ 


\6{fLYE„L 


Cfj — 


\6{fLYE^L 




(3.2) 


m„ here is the multiplicity of momentum states counting the number of states with the same energy 
(mi = 6 ,m 2 = 12, see Ref. [17]). Note that the coefficient c„ has two contributions. One 
stems directly from the 3-pt function, the second one (c^^') has its origin in the denominator of the 
ratio R, i.e. the 2-pt function. 

In figure 2 the coefficients for the two smallest spatial momenta are plotted as functions of 
Mji/Mn for three values of Mj^L. The dependence on these two variables is very mild. Qualitatively 
the same behaviour is found for the coefficients with n > 2 too. The following observations can 
be made: Eirstly, all coefficients are positive. Secondly, is roughly a factor ten smaller than 

hence < 0. And finally, All this implies that the coefficients 
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in the ratio 3.1 are of the same order of magnitude but have opposite sign, bn = bn ^ —Cn- As a 
consequence, the two nucleon-pion-state contributions to the ratio compensate each other to a large 
extent. 

In the following we compare three different ways to obtain an estimate for gA- The first one 
takes g\ K. = t/2) since the excited state contribution is minimal if the insertion time is half 
the sink time. This “midpoint method” is essentially equivalent to the usual plateau method, and 
we refer to it by plateau method 1. The summation method [18] sums R{t,t') over all T between 
source and sink and extracts the axial charge from the term linear in t. This approximation for 
g\ is equivalent to R but with the bn-, -contribution dropped. Finally, one can obtain the axial 
charge also by replacing the 2-pt function in R with the 3-pt function of the vector current (both 
3-pt functions with t' = t/2), and we refer to this as plateau method 2. 

All three methods give estimates for the axial charge that depend on the sink time t. Figure 3 
shows this dependence for a pion mass = ISOMeV and M^L = 4. Plotted is the estimate divided 
by g\, so without the excited state contaminations all three methods give 1. The three curves show 
the result including the nucleon-pion contribution taking into account the first three momentum 
states, i.e. n <3. The yellow region Ifm <t< 1.5fm highlights the range of sink time that seems 
feasible in numerical simulations. The overall conclusion we can draw from this figure is fhaf fhe 
nucleon-pion sfafe confribufion is af fhe few percenf level. For fhe sfandard plafeau mefhod we find 
a 2-3% shiff upwards, while fhe ofher fwo mefhods show a slighfly smaller shifl downwards. The 
resulfs are fairly sfable under changes of fhe pion mass and fhe number of sfafes included in fhe 
confaminafion, see fig. 4. 

Nof obvious from fhese figures is fhe partial compensafion of fhe nucleon-pion-sfafe confribu¬ 
fion menfioned before. For example, for a sink lime t = Ifm in fig. 3 we read off a -i-3% coiTeclion 
for plateau mefhod 1. These 3% are essentially fhe difference (5-2)%, where fhe 5% and 2% stem 
from fhe bn- and c,rConlribulion in R, respectively. If fhe coefficient were approximately equal 
to (and not a factor ten smaller) the nucleon-pion-state contribution would be 5%, so almost a 
factor 2 larger. 

Another consequence of the partial compensation is a fairly flat sink time dependence for 
small pion masses. In the upper left panel of figure 4, for example, fhe fwo plateau methods show 
an almost constant nucleon-pion-state contribution over the yellow sink time region. In contrast, 
the summation method still shows an exponential decay, because the nucleon-pion contribution for 
this method consists of the c„-contribution only and no partial compensation takes place. Whether 
this difference can be resolved in practice is questionable, given that very precise data would be 
needed. Still, the results demonstrate that an increase in sink time not necessarily leads to a smaller 
excited state contribution, at least not for the sink times accessible in practice. 

The results presented here were obtained in a LO calculation. This obviously raises the ques¬ 
tions about higher oder corrections. These are in principle straightforward to compute. However, 
at higher order additional LECs are expected to enter that may be poorly known and the higher 
order corrections are probably less predictive than the LO results discussed here. Another source 
of uncertainty is the mapping of smeared interpolating fields fo poinflike fields in ChPT Naively 
fhis may resulf in an error of 0(M;r/?smear)> which is abouf 30% for a smearing radius of 0.3fm and 
a pion mass less fhan 200MeV. Therefore, if would nof be surprising if fhe LO resulf given here 
has an error of abouf 100%. Lor a more precise esfimafe fhe higher order corrections need fo be 
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Figure 3: Ratio of the estimate and the asymptotic value gx as a function for three different methods to 
obtain the axial charge (see main text). 
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Figure 4: Same as in hgure 3 but with different pion masses (upper two panels) and with different number 
of nucleon-pion-state contributions included (lower two panels). 
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calculated. Still, the sign and the order of magnitude for the nucleon-pion state contribution seem 
to be rather solid results. 

As a final remark we mention that the multi-particle-state contamination with more the one 
pion can be computed in an analogous way. However, the impact of these states is even smaller 
than the nucleon-pion-states considered here. For example, the 3-particle A7r;r-state contribution 
is suppressed by an additional factor [2{fL)'^EjiL]^^ and too small to play a role in practice. 

4. Summary and outlook 

ChPT provides estimates for the nucleon-pion-state contribution in nucleon correlation func¬ 
tions. In case of the nucleon axial charge this contribution is at the few percent level at LO in 
the chiral expansion. Analogous calculations for other observables determined by nucleon 3-point 
functions are straightforward. Particularly interesting are the tensor and scalar nucleon charges, the 
electromagnetic form factors and the quark momentum fraction in the nucleon. 
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